Chern number spins of Mn acceptor magnets in GaAs 
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We determine the effective total spin J of local moments formed from acceptor states bound 
to Mn ions in GaAs by evaluating their magnetic Chern numbers. We find that when individual 
Mn atoms are close to the sample surface, the total spin changes from J = 1 to J = 2, due to 
quenching of the acceptor orbital moment. For Mn pairs in bulk, the total J depends on the pair 
orientation in the GaAs lattice and on the separation between the Mn atoms. We point out that 
Berry curvature variation as a function of local moment orientation can profoundly influence the 
quantum spin dynamics of these magnetic entities. 



State-of-the-art STM techniques have made it possible 
to substitute transition metal impurities for individual 
atoms in semiconductor crystals. [TJ [2] When the impu- 
rity behaves as a dopant, high-resolution STM scanning 
capabilities can then provide detailed information on the 
nature of the bound donor or acceptor states. PUT] Ad- 
vances in spin-polarized STM techniques are now making 
it possible to address the quantum spin dynamics of indi- 
vidual coupled acceptor-impurity systems. [8] These cen- 
ters represent a new class of magnetic entities which we 
refer to as donor magnets or acceptor magnets. They have 
precisely reproducible properties that are intermediate 
in character between those of atomic local moments and 
nanomagnets, have promise for applications in spintron- 
ics and quantum information processing, and act as the 
building blocks of ferromagnetism in semiconductors [9]- 
[TT] that are doped with many transition metal impurities. 

The interpretation of present and future experiments 
requires a theoretical understanding of the quantum dy- 
namics of acceptor and donor magnets. For the specific 
case of individual Mn impurities in bulk GaAs it is known 
that the ground-state total angular momentum of the Mn 
centers is J=l. This value is the result of antiferromag- 
netic coupling between the localized S=5/2 Mn spin and 
the spin (s= 1/2) and orbital moment (1=1) of the ac- 
ceptor hole. The J = 1 character of the ground- state 
effective spin of the Mn embedded atom is supported by 
ESR and infrared spectroscopy experiments. [12} [T3 ] Sev- 
eral questions concerning the magnetic properties of Mn 
dopants in GaAs nevertheless remain unanswered: (i) 
What happens to the total angular momentum J when 
the dopant is close to the symmetry breaking surface 
which provides STM access? (ii) Is there an effective 
"giant spin" describing the low-energy magnetic proper- 
ties of two or more nearby embedded Mn impurities in 
GaAs, and in that case what is its value (iii) Can we 
determine an effective spin Hamiltonian describing the 
quantum dynamics of acceptor magnets? The answers to 
these questions depend on a complex interplay between 
the kinetic exchange that couples Mn and acceptor spins, 
and the variation of the acceptor-level orbital spinor with 



Mn spin orientation which is controlled by spin-orbit in- 
teractions (SOI) and the crystalline environment. 

In this Letter we present a possible answer to these 
questions. The approach we use is similar in spirit to 
ones used to quantize the slow vibronic degrees of free- 
dom in molecular systems. [14 It can be used to quantize 
magnetization dynamics in any theory in which the mag- 
netization direction is initially treated as a classical pa- 
rameter, for example spin-density-functional theory. It 
identifies the effective total spin J of the Mn acceptor 
magnet with a topological Chern number which is the 
average of a Berry curvature electronic functional over 
all possible directions of the Mn acceptor magnetic mo- 
ment. The procedure yields the expected J = 1 value 
for Mn in bulk GaAs. However when a Mn atom is close 
to a symmetry-breaking surface, the circumstance most 
commonly studied in STM experiments, we find that 
J = 2 because the orbital contribution of the acceptor is 
quenched. Mn pairs close to a surface always have J = 4, 
due to strong localization of the acceptor wave-function. 
Surprisingly, for Mn pairs in bulk GaAs we find that the 
total spin can switch between J = 4, J — 3, and J = 2 
depending on the orientation of the pair in the crystal 
and the distance between the two Mn atoms. Our theory 
allows us to extract a quantum spin Hamiltonian for the 
magnetic centers. The spectrum of these Hamiltonians 
can be strongly affected by Berry curvature variation as a 
function of magnetization orientation, which is especially 
strong whenever there is a weakly avoided level crossing 
at the Fermi energy. 

We start by introducing a microscopic tight-binding 
model that captures the salient electronic properties of 
Mn impurities in GaAs. [15j [16] The Hamiltonian reads 

H = i^band + #SO + Jpd ^2 

(i) 

m n [m] 

where i^band contains the Slater-Koster parameters that 
reproduce the band structure of bulk GaAs plus parame- 
ters that account for the As and 4p orbitals of each substi- 
tutional Mn atom; the second term is a one-body atomic 
spin-orbit term. The third term describes an effective an- 
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FIG. 1: (Color online) Electronic structure and Berry curva- 
ture for one Mn impurity in bulk GaAs as a function of the 
polar angle specifying the Mn magnetic moment direction 
relative to a cubic crystal axis, (a) The three acceptor energy 
levels a±, a,2, a% above the valence band edge. The topmost 
level as (red curve) is the only one occupied by a hole; a<i 
(blue curve) is the highest energy occupied orbital, (b) Berry 
curvatures for the individual levels in (a). The dashed lines 
mark the zero for a given curvature plot, (c) Cumulative level 
curvature plus the constant contribution 5/2 from the Mn im- 
purity. The scale for curvature plots is indicated on the right 
axis of panel (c). The blue curve, in which the as contribution 
is not included, is the one that is relevant to the Mn acceptor 
magnet. Adding the curvature of the acceptor level (Z3 to the 
total curvature gives a constant equal to 5/2 (flat red curve). 



tiferromagnetic exchange interaction between a Mn spin 
Sm and its nearest-neighbor As p-spins s n . This is a ki- 
netic pd exchange originating from the hybridization of 
the Mn 3d- levels with the As p- levels. We also include 
a spin-independent Coulombic potential to represent the 
Coulomb potential associated with the Mn ion[15j[T7l[T8]; 
we do not explicitly include the Mn 3d orbit als, and S m 
represents a classical vector of magnitude 5/2. 

The Coulomb potential and level repulsion due to hy- 
bridization with the Mn d-levels together push accep- 
tor states whose spins are aligned with the Mn magnetic 
moments above the valence-band edge.[T5 l \T7 \ H8] Each 
Mn impurity introduces three acceptor levels (p x , p y , p z ) 
which would be degenerate in the absence of SOL SOI's 
not only lift the degeneracy, but also lead to a dependence 
of energies and orbitals on the direction of the Mn mag- 
netic moment. For a neutral Mn impurity only the top 
most of these three states is occupied by a hole. Fig.JIJa) 
shows that this electronic structure is reproduced by the 
model of Eq. [I] implemented numerically for one substi- 
tutional Mn impurity in the middle of a 1200- atom GaAs 
cluster with periodic boundary conditions. For a pair of 
Mn atoms, the lowest energy state of the system is usu- 
ally the one in which the two Mn magnetic moments are 
ferromagnetically aligned. [16] For this configuration, the 
two sets of acceptor states form bonding and antibond- 
ing molecular orbitals. The two topmost empty levels 
are then split by an energy that is related to effective 



Mn-Mn exchange interaction, and varies strongly with 
the pair orientation and the distance between the two 
Mn.[2j[T6] 

This electronic structure sets the stage for the central 
part of our study in which we address the quantum spin 
dynamics of the Mn centers. The starting point to de- 
rive an effective "total spin" Hamiltonian is an approx- 
imate imaginary-time quantum action for the coherent 
spin magnetization direction n[T9] [20] 



S[h] = / dr 



<*[n]|V A *[n]> 



dh 



E[h] 



(2) 



In Eq. ([2]), |\P[n]) is the many-particle ground-state wave 
function obtained by diagonalizing the Hamiltonian ([!]), 
and E[n) is the total energy obtained by summing over 
occupied single-particle states. Here n represents the ori- 
entation of the Mn spin vectors S m , which are assumed 
to be ferromagnetically aligned when the Mn ions are 
not isolated. The first term in ^ is a Berry phase term, 
which for a closed path 7 on the unit sphere is given by 

V = if dh - Vfi^). The line integral can be converted 

into an integral over enclosed area with a gauge- invariant 
integrand known as the Berry curvature [2TJ [22], 



C[h]=iVn X (*|Vfttf). 



(3) 



In our model, C[h] = 5^° cc Ci[n], the sum of the Berry 
curvatures Ci[h] of all occupied single-particle levels i = 
1,2, ...,occ. In the absence of SOI, all C{ = Ci[n] ■ h 
are constant and equal to the spin projection of each 
level, ±1/2. In this case C • h = ^\ ~ $i where 

S is the total spin of the system. The effect of SOI is 
twofold. First, there is now an orbital contribution to 
the curvature. Second, Ci[h] varies with n, as shown in 
Fig. [TJb). 

The average of the Berry curvature over the unit sphere 
of all possible directions, 



J 



Js 2 



C[f 



~n dA 



(4) 



is a topological invariant, known as first Chern 
number. [23] Although C[n] can fluctuate strongly, J is 
always a half-integer; its value can change only if the 
system suffers a level crossing at the Fermi level. [T9] [20] 
It is useful to introduce single-level Chern numbers ji = 
^ J S 2 Ci[h] • n dA, so that J = ^° cc When SOI's are 
included ji can be different from ±1/2. The contribu- 
tion from the spin-polarized Mn d-orbitals that are not 
explicitly included in our calculation are accounted for 
by adding 5/2 to the total Berry curvature for each Mn 
ion in the system. As we explain below, the total Chern 
number J plays the role of the effective quantum spin of 
the Mn acceptor magnetic. 

We first discuss the case of one Mn impurity in GaAs 
and confirm that our approach yields sensible results. 
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Indeed, as shown in Table I, when the Mn is located 
in the bulk, solving our model numerically for a 1200- 
atom cluster gives J = 1, as expected. To understand 
this result, we observe that the total Chern number for 
a full valence band vanishes. The total Chern number of 
the acceptor magnet is therefore J = 5/2 — j, where 
j is the Chern number of the topmost acceptor state 
which is occupied by a hole. As shown in Table I, we 
find j = 3/2 = 1/2 + 1. The orbital contribution to j, 
which is absent in any theory that does not include spin- 
orbit interactions appears in our approach because the 
orbital content of the topmost acceptor state varies with 
moment orientation. The orbital moment is locked to 
the spin-moment by strong SOI's. The result is different 
however, when the Mn impurity is located on the (110) 
surface, [24] a case frequently considered in STM exper- 
iments. For this case we find that the orbital content 
of the acceptor wavefunction does not vary substantially 
with magnetization direction. The Chern number of the 
acceptor level only has a j = 1/2 spin contribution so 
that J = 5/2 — 1/2 = 2. Evidently symmetry breaking 
on the (110) surface, with only three nearest-neighbor As 
atoms instead of four, creates a local environment for a 
Mn whose symmetry is lower than the tetragonal sym- 
metry seen by the impurity in bulk GaAs. As a result the 
acceptor orbital moment is completely quenched. (A di- 
rect calculation shows that (as\L\as) =0.) Interestingly, 
our calculations show that j remains equal to 1/2 also 
when the Mn is located on one of the immediate sub- 
surfaces below the top (110) surface. We expect that 
j should eventually switch to the bulk value when the 
Mn is located deeply below the surface, but this does not 
happen for the film thicknesses in the present simulation. 

We can now investigate magnetic clusters with two 
Mn atoms, where the resulting total spin is less intu- 
itive. Table II shows the Chern numbers for several 
ferromagnetic Mn pairs, whose orientation (Iran) in the 
crystal and Mn separation d is described in Fig. [5|c). 
When the pair is positioned on the (110) surface or in 
one of its nearby sub-layers, we find that the individ- 
ual Chern numbers of the two empty acceptor states are 
always j\ = 22 — 1/2, yielding a total Chern number 
7 = 2*5/2-2*1/2 = 4. The situation for a Mn pair in 
bulk GaAs is more complex. We can see that, while for 
several pairs (e.g. all the (110) pairs) J = 4, like for the 
surface case, other pairs have J = 3, and ferromagnet- 
ically coupled remote spins should have J = 2. In this 
second case one of the individual acceptor Chern num- 
bers, ji or j'2, is equal to 3/2, and in the latter case both 
have Chern number 3/2. Clearly the presence of a second 
Mn affects the orbital magnetic properties of the other, 
in a way that depends both on pair orientation and Mn 
separation. The outcome for the pair is not easily pre- 
dictable. For example, while the (100) pair switches from 
J — 3 at the shortest Mn separation d = a to J = 4 at 
larger separations, the (211) pair behaves exactly in the 
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FIG. 2: (Color online) Calculated Chern numbers for Mn 
impurities in bulk GaAs or in the (110) surface and nearby 
sub-surfaces, (a) Results for one Mn impurity. J is the total 
Chern number and j the Chern number of the acceptor level, 
(b) Results for various Mn pairs (lmn)d separated by distance 
d in units of the lattice constant a. j\ and are the two ac- 
ceptor Chern numbers, (c) Orientation (Imn) and separation 
of different Mn pairs on the (110) surface of GaAs. 



opposite way. (a is the lattice constant.) 

In the remaining part of the paper we will extract a 
quantum spin Hamiltonian describing the dynamics as- 
sociated with the moment orientation of acceptor mag- 
nets. We return to the action given by Eq. [2] and 
perform a change of variables p~9j [20] from n (#,(/>) to 
h f (6 f ,(j)') that transforms the Berry curvature C[h] to a 
constant C'[n'] = J. This change of variables rescales 
the local curvature metric such that C (#, 0) s'm(6)d0d(j) = 
J sm(0)d0' ' d<j)' . The real-time action for a path becomes 



/' 

Jo 



dt' 



i A, 



E{h[ti 



(*')]} 



(5) 



where Aj = J<p' (1 - cos6>') /sin6>'. Eq. |] is the quan- 
tum action for an effective "total spin" quantum number 
J[22 . The second term in the integrand is the semi- 
classical Hamiltonian of the system, which is given by 



E{n') = (E{n[h']} = (J,h'\H\J,n') 



(6) 



where % is the quantum spin Hamiltonian and |J, n') is 
a spin-J coherent state parametrized by the unit vector 
n'(6' , (j)'). The function E(n') is is the anisotropy energy 
transformed so that it also captures Berry curvature vari- 
ation. The quantum Hamiltonian Iri is constructed by 
first expanding E(n') in spherical harmonics E(n') — 
Y?tUYLL=-n? Y r^'). We then use a formula[20] 
which relates the spherical harmonics expansion coeffi- 
cients 7^ to the matrix elements = (J,m\i-L\J,m') 
of the quantum spin Hamiltonian: 



H 



2J X 

(j) -(_i r '-' ^ v; ^ 



A=0 fi=-\ 



2A+ 1 



47T 



J 

rn 



J A 
-m' [i 



J J X 
J -J 



(7) 



4 



ft (211)d =2 .4g C{9 (j)) EijmeV] 

3 ^ ' ' ' 'l^SHPn'^g 30 ~° ' ' ' ' — 

%^^4 5 6 f™ 1 °2 V 6 V 

(a) 9 (b) (EilS^m 

n( 211) d=2Aa (H^)[meV ] i2 z ft (211), =2 , 4a (^ (3) ) [meV ? ? 

r : ® <fc ! > 1 *[ | 

1 2 3 4 5 6a 1 2 3 4 5 6a 
(C) ^ (d) 9 

FIG. 3: (Color online) Effects of the Berry curvature on 
the quantum spin properties of the ferromagnetic Mn pair 
(211)d = 2.4a in bulk GaAs. (a) Total Berry curvature land- 
scape as a function of the pair magnetic moment direction. 
The two large dips at ~ 0, 7r, relative to the [001] axis signal 
avoided crossings at the Fermi level, (b) Spectrum of the ef- 
fective quantum spin Hamiltonian with Chern number J — 3 
vs. the expectation value of The filled circles are the 
spectrum of 7i 3 obtained by quantizing the classical magnetic 
anisotropy landscape. The empty circles are the spectrum of 
Ti 3 , which includes Berry curvature effects, (c) and (d) are 
the magnetic anisotropy landscape for the cases when Berry 
curvature effects are excluded or included respectively. 



tal spin of Mn impurities in GaAs is a topological Chern 
number that includes both the contribution of the Mn 
spins and the spin and orbital moments of the acceptor 
states. The effective spin depends sensitively on the en- 
vironment around the Mn impurities and it is strongly 
affected by the presence of symmetry breaking surfaces 
and the geometry of the magnetic clusters in the GaAs 
lattice. The quantum dynamics of the effective spin is 
qualitatively modified by Berry phase corrections caused 
by electronic degeneracies at some direction of the Mn 
magnetic moments. This has implications on the stabil- 
ity of the quasi-classical magnetization of the magnetic 
center around minima in the anisotropy energy. 
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where the quantities in parenthesis are Wigner 3J sym- 
bols. Once the Hamiltonian matrix has been obtained, 
it can be decomposed and rewritten as a combination of 
spin operators. [20] 

As an example of this procedure, we show in Fig. [^re- 
sults for the Mn pair (211)d=i.4 a m bulk GaAs. In panel 
(a) we plot the Berry curvature functional C[n] = C(0, 
The coordinate system used for these plots has = par- 
allel to the [001] axis, (0 = tt/2,0 = 0) parallel to [100], 
and (0 = tt/2, <j> = tt/2) parallel to [010] (see Fig. [2). The 
calculated Chern number for this pair is J = 3. The large 
dips of C(0, qb) below this value for « 0, 7r signal the oc- 
currence of narrowly avoided level crossings at the Fermi 
level, for two time-reversed directions. According to our 
theory, we expect that in this case Berry phase vari- 
ations strongly influence the spectrum of the quantum 
spin Hamiltonian. This is indeed the case, as shown in 
Fig. |3|b), where we plot the spectrum {£^}i=i,...,2j+i of 
the Hamiltonian obtained when Berry phase corrections 
are either included (H 3 , empty circles) or absent (H 3 , 
filled circles), versus (Ei\S^\Ei) . (S z is the z-component 
of the effective spin.) The difference in the two spec- 
tra has implications for the classical magnetic anisotropy 
landscape. In Fig(3^c) and (d) we plot (J, n|H 3 |J, n) and 
( J, n\H 3 \ J, n) respectively, as a function of and <ft. Al- 
though the anisotropy minima are present in both cases, 
the barrier that separates them is considerably reduced 
by Berry phase corrections. The dips in the curvature at 
level crossings increase the quantum tunneling rates of 
the magnetization between the two minima. 

In conclusion, we have proposed that the effective to- 
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